Confinement of electrons in layered metals 
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We analyze the out of plane hopping in models of layered systems where the in-plane properties 
deviate from Landau's theory of a Fermi liquid. We show that the hopping term acquires a non 
trivial energy dependence, due to the coupling to in plane excitations, and can be either relevant or 
irrelevant at low energies or temperatures. The latter is always the case if the Fermi level lies close 
to a saddle point in the dispersion relation. 



Layered materials have been object of intensive study 
since they present important physics. Unusual properties 
are derived from the anisotropy and periodicity along 
the axis perpendicular to the planes Among the 

most studied layered materials are the high-temperature 
cuprate superconductors. These compounds present a 
strong anisotropy and are treated as two-dimensional sys- 
tems in many approaches. In the normal state the trans- 
port properties within the Cu02 planes are very differ- 
ent from those along the c-axis: electron motion in the 
c-direction is incoherent in contrast with the metallic be- 
havior of the in-plane electrons as probed by the dif- 
ferent pc and pab resistivities and their different depen- 
dence with temperature [^,|3| . Optical conductivity mea- 
surements confirm the anomalous c-axis properties [Q. 
The relevance of the nature of the conductance in the 
direction perpendicular to the Cu02 planes to the na- 
ture of the superconducting phase has been remarked on 
both theoretical and experimental |^ grounds. The 
anomalous nature of the out of plane properties in the 
cuprates, and in analogy with the one dimensional Lut- 
tinger liquid, has led to the proposal of the failure of the 
conventional Fermi- liquid theory in these compounds . 
An alternative explanation of the emergence of incoher- 
ent behavior in the out of plane direction has been pro- 
posed in terms of the coupling of the interlayer electronic 
motion to charge excitations of the system [Q. This ap- 
proach implicitly assumes that electron electron interac- 
tions modify the in-plane electron propagators in a non 
trivial way, at least at distances shorter than the elas- 
tic mean free path, I. The strong angular dependence 
of both the scattering rate and the interplane tunneling 
element can also lead to the observed anisotropics 

Graphite, another layered material, presents a in- 
traplane hopping much larger than the interplane hy- 
bridization. The unconventional transport properties of 
graphite such as the linear increase with energy of the in- 
verse lifetime |Q (see also ) , suggest deviations from 
the conventional Fermi liquid behavior, which could be 
due to strong Coulomb interactions unscreened because 
of the lack of states at the Fermi level 



In the following, we extend the methods used in 0| to 
the clean limit, Z — > oo, and show that, for certain models 
of correlated electrons the interplane hopping be- 

tween extended states can be a relevant or an irrelevant 
variable, in the RG sense, depending on the strength of 
the couphng constant. The scheme used here is based on 
the renormalization group analysis as applied to models 
of interacting electrons . For the problem of inter- 

chain hopping between Luttinger liquids, it recovers well 
known results ]T^-p^. 

We first present the method of calculation, and show 
the results it gives in one dimension. Then, we apply it 
to two dimensional models which show deviations from 
FL behavior. The main physical consequences of our 
calculation are discussed in the last section. 

The method of calculation. In the presence of electron- 
electron interactions, tunneling processes are modified 
by inelastic scattering events. If the excitations of the 
systems (electron-hole pairs, plasmons) are modelled as 
bosonic modes, one can write an effective electron-boson 
hamiltonian of the type: 



\ ^kblbk + fffc,»4c,(&^ -I- bk) 

(1) 



where the ty describe the electronic hopping processes. 
The electron-boson interaction leads to a Franck Condon 
factor which reduces the effective tunneling rate. The 
electron propagators acquire an anomalous time, or en- 
ergy, dependence: 
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where: 
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This expression can be generalized, taking into account 
the spatial structure of the coupling to: 
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where J7 is the region of overlap of the wavefunctions 
\l/i(f) and 4'j(r). This expression, which can be seen as 
the exponential of the leading frequency dependent self- 
energy correction to the electron propagator, has been 
extensively used in studies of tunneling in zero dimen- 
sional systems (single electron transistors) which show 
Coulomb blockade , one dimensional conductors , 
and disordered systems in arbitrary dimensions pl| ] . 

The effective interaction can, in turn, be written in 
terms of the response function as: 



(5) 



In a Fermi liquid, we have Ve//(k,Ci;) w Q!(k)|w| for 
uj ^ Ep, where Ep is the Fermi energy, so that: 

lim (*t(t)^(t'))^ 
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The parameter a gives the correction to the scaling prop- 
erties of the Green's functions. It is easy to show that, in 
an isotropic Fermi liquid in D dimensions, lim/^oo a oc 
where I is the linear dimension of the (local- 
ized) electronic wavefunctions 5'(f). This result is due 
to the dependence on k of the response functions, as 
Imx(k,w) ^ |w|/(fc^~^|k|). Thus, for D > 1, we recover 
coherent tunneling in the limit of delocalized wavefunc- 
tions. 

In one dimension, one can use the non interacting ex- 
pression for Imxo(k, w), to obtain: 
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t-t' <^ l/vp 



^ log[vp{t-t')/l] t~t'^l/vp 



where we have assumed a short range interaction, U. 
Hence, the Green's functions have a non trivial power 
dependence on time, even in the / — > oo limit, in agree- 
ment with well known results for Luttinger liquids [po[ . 
In order to obtain the energy dependence of the effec- 
tive tunneling between k states near the Fermi surface, 
one needs to perform an additional integration over dr. 
In general, near a scale invariant fixed point, uj oc |k|^, 
and for a ID conductor one knows that z = 1. Hence, 



ImG(tJ, kp) (X to' 



-z+a 



The flow of the hopping 



terms under a Renormalization Group scaling of the cut- 
off is ITT 
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—a localized hopping 
1 — a extended hopping 



(8) 



where t denotes a hopping term, between localized or 
extended states. In the latter case, the hopping becomes 
an irrelevant variable for a > 1. 

Graphene planes. The simplest two dimensional model 
for interacting electrons where it can be rigourously 
shown that the couplings acquire logarithmic corrections 
in perturbation theory is a system of Dirac fermions 
(cfc — vf\^\), with Coulomb, l/|f — r'|, interaction. This 
model can be used to describe isolated graphene planes 
p2[ , and can help to understand the anomalous behavior 
of graphite observed in recent experiments jlO[ ^ . 

In order to apply the procedure outlined in the previ- 
ous section, one needs the Fourier transform of the inter- 
action, e^/(eo|k|), where e is the electronic charge, and 
eo is the dielectric constant, and the susceptibility of the 
electron gas. For a single graphene plane, this quantity 
is: 



Xo(k,cj) 



1 



(9) 



These expressions need to be inserted in equations (|^) 
and (^). Alternatively, we can use the RPA, and include 
the effects of interplane screening: 



XflPA(k,cj) = 



ih(|k|d) 



[cosh(|k|d) + |^sinh(|k|<i)xo(k.t^)] 



(10) 



where d is the interplane spacing. The imaginary part, 
Imxo(k,ct;), is different from zero if w > WF|k|. 

For simplicity, we consider the expression in eq.(|^), 
as it allows us to obtain analytical expressions. We cut 
off the spatial integrals at a scale, I, of the order of the 
electronic wavefunctions involved in the tunneling. We 
obtain an expression similar to that in eq.(|^) except that 
the prefactor {U/Ep)"^ is replaced by e'^/{eoVp)'^. Thus, 
the propagators acquire an anomalous dimension. As in 
ID, the value of the exponent z which relates length and 
time scales is z = 1. The scaling of the hoppings now 
are: 
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The extra constant in the first equation with respect to 
eq. (^) reflects the vanishing of the density of states at 
the Fermi level for two dimensional electrons with a Dirac 
dispersion relation. 

In graphite, the dimensionless coupling constant, 
/vp, is of order unity. Under renormalization, it flows 
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towards zero [Q. Thus, interplane tunneling is a rele- 
vant variable, although, in a dirty system with a finite 
mean free path, interplane tunneling can become irrele- 
vant n. 

Saddle point in the density of states. The Fermi sur- 
face of most hole-doped cuprates is close to a Van Hove 
singularity . The possible relevance of this fact to the 
superconducting transition as well as to the anomalous 
behavior of the normal state was put forward in the early 
times of the cuprates and gave rise to the socalled Van 
Hove scenario p^ . 

We shall apply the mechanism described in section II 
to study the interlayer hopping of two electronic systems 
described by the t— t' Hubbard model and whose Fermi 
surface lies close to a Van Hove singularity. 

The t-t '-Hubbard model has the dispersion relation 



e(k) = — 2t [cos{kxa) + cos(fcj,a)] 
— 2t' cos(/ca;a) cos(fc.ya) , 



(12) 



This dispersion relation has two inequivalent saddle 
points at A (7r,0) and B (0,7r). The Van Hove model 
in its simplest formulation is obtained by assuming that 
for fillings such that the Fermi line lies close to the sin- 
gularities, the majority of states participating in the in- 
teractions will come from regions in the vicinity of the 
saddle points |2|] . Taylor expanding eq. ( p^ ) around the 
two points gives the effective relation 



saM'^) ~ T(t T 2t')kl ± (t ± 2t')fc^ 
which leads to the effective hamiltonian: 



(13) 
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where Ag is a high energy cutoff which sets the limit of 
validity of the effective description. The particle-hole 
susceptibility has been computed in |13[|: 



Imx(k,Lj) = 



1 



47r£(k) 



e(k) - u;-e(k) 



(15) 



where e(k) is the dispersion relation (|l3|). 

The long time dependence of the Green's function is 
determined by the low energy behavior of x- 



lim Im x(k, uj) 

UJIU — >0 



e(k) 



Inserting this expression in eqs.(y) and (g), we can see 
that, irrespective of the details of the interaction, in the 
presence of a Van Hove singularity the exponent a in the 
time dependence of the Green's function goes as: 

lim a (X \og{l) , 



where I, as before, is the length scale which characterizes 
the wavefunction of the tunneling electron, and 17 cx P is 
the size of the integration region in cq.(^). 

The details of the anomalous dimension of the prop- 
agator in this case depend on the nature of the interac- 
tions which determine Ve//(k, cj) — V^^(k)Imx(k, w) . To 
make contact with previous works 0| , we have computed 
the parameter a for different possible interactions: 

a) Unscreened Coulomb potential, VCk) — Due 

Co |k 

to the highly singular interaction, the uj dependence of 
the effective potential is not linear and a is not well de- 
fined. The effective potential eq. (||) is computed to be: 



VeffiKuj) 







1 -l-log 


(wo)_ 







where ojq = \Ep — Evh\ is a low-energy cutoff, required 
to avoid divergences in the integrals. The Fermi energy 
is Ep and the position of the saddle point is at Evh- 

b) Screened interaction of the Hubbard type, V^(k) = 
Ua^, where a is the lattice unit. 
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7rt2 (27r)2 
where A is a high-energy cutoff, and 
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where kx,ky ~ ^/oJo■ 

c) Thomas-Fermi screened potential. 



(2^e2)2 



kprp) 



2c k jpip 
eo7rt2 



log'ikFT/k) 



where kpx is the Fermi- Thomas wavevector, and, as be- 
fore, k ~ ^yuJo■ This is an intermediate situation between 
the two previous cases. 

Finally, in the case of a layered system, we can use 
the RPA summation, including interlayer interactions, 
eq.(|l^). In this case, z — 2, and the scaling equations for 
the hoppings are the same as in eq.(|^). 

A numerical computation provides the effective poten- 
tial shown in fig. |l|. 




FIG. 1. Effective potential eq. (|^ as a function of the 
energy for fixed k for a system of Van Hove layers coupled by 
Coulomb interaction . 
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Fig. |l|] compares well with the experimental plots 
of the loss function given in what reveals that the 
Van Hove model is also compatible with transport exper- 
iments. The results for the parameter a extracted from 
the numerical computation are in qualitative agreement 
with the analytical expressions given in cases b) and c) 
above. 

In all cases a diverges as Ep Evh- Thus, interlayer 
hopping is an irrelevant variable, and scales towards zero 
as the temperature or frequency is decreased. This addi- 
tional logarithmic dependence can be seen as a manifesta- 
tion of the log^ divergences which arise in the treatment 
of this model jl^. Note that, as in the graphene case, the 
coupling constants are also energy dependent, and grow 
at low energies, suppressing even further the interlayer 
tunneling. 

Conclusions. We have discussed the supression of in- 
terlayer tunneling by inelastic processes in two dimen- 
sional systems in the clean limit. Our results suggest 
that, when perturbation theory for the in-plane interac- 
tions leads to logarithmic divergences, the out of plane 
tunneling acquires a non trivial energy dependence. This 
anomalous scaling of the interlayer hopping can make it 
irrelevant, at low energies, if the in-plane interactions are 
sufficiently strong. This is always the case if the Fermi 
level of the interacting electrons lies at a van Hove singu- 
larity (note that the Fermi level can, in certain circum- 
stances, be pinned to the singularity [^). Thus, we have 
shown that insulating behavior in the out of plane direc- 
tion is not incompatible with gapless in-plane properties. 
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